Novel evaluation of the two-pion contribution to the nucleon 

isovector form factors 



M.A. Belushkiifl and H.-W. Hammei0 

Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie), 
Universitdt Bonn, Nufiallee 14-16, D-53115 Bonn, Germany 

Ulf-G. MeifineH 

Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie), 
Universitdt Bonn, Nufiallee 14-16, 
D-53115 Bonn, Germany 
and 

Institut fur Kernphysik (Theorie), 
Forschungszentrum Jiilich, 
D-52425 Jiilich, Germany 

(Dated: February 2, 2008) 

Abstract 

We calculate the two-pion continuum contribution to the nucleon isovector spectral functions 
drawing upon the new high statistics measurements of the pion form factor by the CMD-2, KLOE, 
and SND collaborations. The general structure of the spectral functions remains unchanged, but 
the magnitude increases by about 10%. Using the updated spectral functions, we calculate the 
contribution of the two-pion continuum to the nucleon isovector form factors and radii. We compare 
the isovector radii with simple p-pole models and illustrate their strong underestimation in such 
approaches. Moreover, we give a convenient parametrization of the result for use in future form 
factor analyses. 
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I. INTRODUCTION 



The electromagnetic form factors of the nucleoli offer a unique window on strong in- 
teraction dynamics over a wide range of momentum transfers 0, 0|. At small momentum 
transfers, one is only sensitive to the gross properties of the nucleon like the charge and 
magnetic moment, while at high momentum transfers one will be able to resolve aspects of 
the quark substructure of the nucleon as described by QCD. The form factors encode both 
perturbative and non-perturbative aspects of QCD such as the nucleon radii, vector meson 
coupling constants, and the asymptotic behavior at large momentum transfer. Moreover, 
their detailed understanding is important for a wide variety of experiments ranging from 
the strange vector form factors of the nucleon [3| to Lamb shift measurements in atomic 
hydrogen jj. 

It has been known for a long time that the pion plays an important role in the long- 
range structure of the nucleon [5(. This connection was made more precise using dispersion 
theory in the 1950's 0,0]. Subsequently, Frazer and Fulco have written down partial wave 
dispersion relations that relate the nucleon electromagnetic structure to pion-nucleon (ttN) 
scattering and predicted the existence of the p resonance jH || . Despite of this success, the 
central role of the 2tt continuum in the isovector spectral function has often been ignored in 
vector-meson dominance analyses of the electromagnetic form factors of the nucleon where 
the p was approximated by a simple pole. In 1975, Holder and Pietarinen pointed out that 
this omission leads to a gross underestimate of the isovector radii of the nucleon [10]. Holder 
et al. first perfomed a consistent dispersion analysis of the electromagnetic form factors of 
the nucleon including the 2tt continuum derived from the pion form factor and nN- 
scattering data [l2l |. In the mid-nineties, this analysis has been updated and extended 
to include data in the time-like region |lj]. Recently, the new precise data for the neutron 
electric form factor have been included as well 2 • 



Using chiral perturbation theory (CHPT), the long-range pionic structure of the nucleon 
can be connected to the Goldstone boson dynamics of QCD ^6|. The nonresonant part 
of the 2ir continuum is in excellent agreement with thephenomenological analysis and 
the p- meson contribution can be included as well 0, 0, It is well known that vector 
mesons play an important role in the electromagnetic structure of the nucleon, see e.g. 
IS Q Q UBS, and the remaining contributions to the spectral function have 



usually been approximated by vector meson resonances. 

A new twist to this picture was recently given in Ref. [26|, where the form factor data 
were interpreted based on a phenomenological fit with an ansatz for the pion cloud using the 
old idea that the proton can be thought of as virtual neutron-positively charged pion pair. 
A very long-range contribution to the charge distribution in the Breit frame extending out 
to about 2 fm was found and attributed to the pion cloud. While naively the pion Compton 
wave length is of this size, these findings are indeed surprising if compared with the "pion 
cloud" contribution due to the 2n contribution for the isovector spectral functions of the 
nucleon form factors, which can be obtained from unitarity or chiral perturbation theory. 
As was shown in Ref. (2?}, these latter contributions to the long-range part of the nucleon 
structure are much more confined in coordinate space and agree well with earlier (but less 
systematic) calculations based on chiral soliton models, see e.g. |28[. Therefore it remains 
to be shown how to reconcile the findings of Ref. |26[, based on a global fit to all nucleon 
form factors, with the results of dispersion analysis and chiral perturbation theory. 

The CMD-2 j2^|, KLOE [3(j, and SND j3]J collaborations have recently remeasured the 



2 



pion form factor with high statistics and found significant deviations from earlier measure- 
ments at the p-resonance peak. Moreover, the three measurements also show discrepancies 
among each other. In light of the importance of the 2ir continuum for the electromagnetic 
structure of the nucleon, it appears worthwhile to recalculate the 27r-continuum contribu- 
tion to the isovector spectral function of the nucleon using the new high statistics data and 
estimate the errors that arise from the discrepancy among the experiments. Moreover, a 
new analysis will help to better understand the nature and range of the pion cloud of the 
nucleon. 



II. NUCLEON FORM FACTORS 

First, we collect some basic definitions. The nucleon electromagnetic (EM) form factors 
are defined by the nucleon matrix element of the electromagnetic current, 



W)\J™\N(p))=u(p') 



^ Fl (t) + ^a^(p'-p) v F,(t) 
2m 



u(p) , (1) 



with t = (p' — p) 2 = q 2 < the invariant momentum transfer squared, and m the nucleon 
mass. Fi(t) and F2(t) are the Dirac and the Pauli form factors, respectively. They are nor- 
malized at t — to the charge (Fi) and anomalous magnetic moment (F 2 ). It is convenient 
to work in the isospin basis and to decompose the nucleon form factors into isoscalar (S) 
and isovector (V) parts, 

F?=l(F?+Fn, pr = 1 -(F[-Fn, 1 = 1,2, (2) 

where p (n) denotes the proton (neutron). The experimental data are usually given for the 
Sachs form factors, which are linear combinations of F\ and F 2 

G I E (t) = F((t)-rFi(t), G I M (t) = F((t) + F 2 I (t), I — S,V , (3) 

where r = — t/(4m 2 ). In the Breit frame, where the energy transfer of the virtual photon 
vanishes, Ge and Gm may be interpreted as the Fourier transforms of the charge and 
magnetization distributions in coordinate space, respectively. 

The analysis of the nucleon electromagnetic form factors proceeds most directly through 
the spectral representation given by 1 

^ i r^m*, i = 1 , 2 , i=SiVi (4) 

71 Ifl t — t 

where the corresponding thresholds are given by t$ = (3M W ) 2 and t% = (2MJ 2 , respectively 
and M n is the pion mass. The imaginary part entering Eq. (|3j) can be obtained from a 



1 We work here with unsubtracted dispersion relations. Since the normalizations of all the form factors are 
known, one could also work with once-subtracted dispersion relations. For the topic studied here, this is 
of no relevance. 
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spectral decomposition For this purpose it is convenient to consider the EM current 
matrix element in the time-like region, 



J, = (N(p)N(p)\J^\0)=u(p) 



F 1 (t) 1 , + ^a, u (p + pyF 2 (t) 



v(p) , (5) 



where p , p are the momenta of the nucleon-antinucleon pair created by the current J^ m . The 
four-momentum transfer in the time-like region is t = {p + p) 2 > 0. Using the LSZ reduction 
formalism, the imaginary part of the form factors is obtained by inserting a complete set of 
intermediate states as p| 

ImJ, = |(27r) 3 / 2 Ar^(p|J JV (0)|A)(A|J7|0)^)5 4 (p + P-PA), (6) 

A 

where M is a nucleon spinor normalization factor, Z is the nucleon wave function renor- 
malization, and Jn{x) = J^(x)7o with Jn{x) a nucleon source. The states |A) are asymp- 
totic states of momentum p\. They carry the same quantum numbers as the current J^ m : 
I G (J PC ) = 0~(1 ) for the isoscalar component and I G (J PC ) = 1 + (1 ) for the isovector 
component of J^ m . Furthermore, they have no net baryon number. 

For the isoscalar current the lowest mass states are: 3tt, 5tt, . . .; for the isovector part they 
are: 2tt, Att, .... Because of G-parity, states with an odd number of pions only contribute 
to the isoscalar part, while states with an even number contribute to the isovector part. 
Associated with each intermediate state is a cut starting at the corresponding threshold in 
to and running to infinity. As a consequence, the spectral function 1m F(t) is different from 
zero along the cut from t to oo with t = 4 (9) M 2 for the isovector (isoscalar) case. Using 
Eqs. (|5I6|) . the spectral functions for the form factors can in principle be constructed from 
experimental data. In practice, this program can onl y b e carried out for the lowest-mass 
two-particle intermediate states {2n and KK) [lil. I32T.I33T]. 

The longest-range (and therefore most important at low momentum transfer) pion cloud 
contribution comes from the 2-k intermediate state in the isovector form factors. This con- 
tribution was first constructed from the pion form factor and ttN scattering data by Holder 
and Pietarinen llS 



III. TWO PION CONTINUUM 

In this paper, we re-evaluate the 2-ir contribution in a model-independent way using the 
latest experimental data for the pion form factor from CMD-2 [2^, KLOE j^jj, and SND |3l| . 
We follow Ref. 0] and express the 2n contribution to the the isovector spectral functions 
in terms of the pion charge form factor F n (t) and the P-wave irn — > NN amplitudes fl(t). 
The 2ir continuum is expected to be the dominant contribution to the isovector spectral 
function from threshold up to masses of about 1 GeV 34]. Here, we use the expressions 

ImG^(t) = -^-F„(t)* fl(t) , 
m\Jt 

ImG^(t) = -^F^tyf^t), (7) 



where q t = y/t/4 - M%. The imaginary parts of the Dirac and Pauli Form factors can 
be obtained using Eq. (jSJ). The P-wave nn — ► NN amplitudes f±(t) are tabulated in 
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t [GeV 2 ] 

FIG. 1: The pion electromagnetic form factor in the time-like region as a function of the momentum 
transfer t. The diamonds, squares, and circles show the high statistics data from the CMD-2 |2sj |. 
KLOE [3(|, and SND [ijj collaborations, respectively. The dashed, solid, and dash-dotted lines 
are our model parametrizations. The inset shows the discrepancy in the resonance region in more 
detail. 



Ref. [34j. (See also Ref. [35[ for an unpublished update that is consistent with Ref. |34(.) 
We stress that the representation of Eq. (J7J) gives the exact isovector spectral functions for 
4M 2 . < t < 16M 2 , but in practice holds up to t ~ 50M 2 . Since the contributions from An 
and higher intermediate states is small up to t ~ 50M 2 , F n (t) and the f±(t) share the same 
phase in this region and the two quantities can be replaced by their absolute values. 2 

The updated pion form factor is shown in Fig. ^ The diamonds, squares, and circles 
show the high statistics data from the CMD-2 [29], KLOE j3fj, and SND [31| collaborations, 
respectively. The dashed, solid, and dash-dotted lines are our model parametrizations which 
are of the Gounaris-Sakurai type r! 2^1- The form factor shows a pronounced p-u mixing 



in the vicinity of the p-peak. There are discrepancies between the three experimental data 
sets for the pion form factor j3jj . The discrepancies in the p-resonance region are shown in 
more detail in the inset of Fig. ^ Since we are not in the position to settle this experimental 
problem, we will take the three data sets at face value. We will evaluate the 27T continuum 
given by Eq. (|7J) for all three sets and estimate the errors from the discrepancy between the 
sets. 



Using the new high statistics pion form factor data [29J, |3CJ, |31| and the amplitudes 



2 We note that representation of Eq. Q) is most useful for our purpose. The manifestly real functions 
J±(t) = /_J_ (t) /i 7 ^ (0 also tabulated in Ref. |2J| contain assumptions about the pion form factor which 
leads to inconsistencies when used together with the updated F v (t). 
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FIG. 2: The 2tt spectral function using the new high statistics data for the pion form factor 
[31( 1 . The spectral functions weighted by 1/t 2 are shown for Ge (solid line) and Gm 
(dash-dotted line). The previous results by Hohler et al. [3J| (without p-uj mixing) are shown for 
comparison by the gray/green lines. 



f±(t) tabulated in Ref. |34j . we obtain the spectral functions shown in Fig. El We show 
the spectral functions weighted by 1/t 2 Ge (solid line) and Gm (dash-dotted line). The 
previous results by Hohler et al. (without p-u mixing) are given for comparison by the 
gray/green lines. The general structure of the two evaluations is the same, but there is a 
difference in magnitude of about 10%. The difference between the three data sets for the 
pion form factor is very small and indicated by the line thickness. The difference in the 
form factors is largest in the p-peak region (cf. Fig. |TJ), but this region is suppressed by the 
7T7T — > NN amplitudes f±(t) which show a strong fall-off as t increases. 

The spectral functions have two distinct features. First, as already pointed out in |8|, they 
contain the important contribution of the p-meson with its peak at t ~ 30M 2 . Second, on the 
left shoulder of the p, the isovector spectral functions display a very pronounced enhancement 
close to the two-pion threshold. This is due to the logarithmic singularity on the second 
Riemann sheet located at t c = AM 2 - M*/m 2 = 3.98M 2 , very close to the threshold. This 
pole comes from the projection of the nucleon Born graphs, or in modern language, from 
the triangle diagram. If one were to neglect this im por tant unitarity correction, one would 
severely underestimate the nucleon isovector radii |l0j. In fact, precisely the same effect 
is obtained at leading one-loop accuracy in relativistic chiral perturbation theory 
This topic was further elaborated on in the framework of heavy baryon CHPT [17, 20] and 
in a covariant calculation based on infrared regularization |l8|. Thus, the most important 
27r contribution to the nucleon form factors can be determined by using either unitarity or 
CHPT (in the latter case, of course, the p contribution is not included). 

The importance of the correct 27r-continuum contribution for the nucleon isovector radii 
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T > r 

Kei. 


(r )[ |im \ 




this work 


0.32.. .0.33 


1.77.. .1.80 


Ref. [20] 


0.19 


1.10 


Ref. f38] (*) 


0.19 


0.96 


Ref. [38] (u) 


0.27 


1.38 



TABLE I: Comparison of the 27r-continuum contribution to the Dirac and Pauli isovector radii 
of the nucleon (first line) to three different p-pole models from Refs. Big. The given range 
indicates the error introduced by the different pion form factor data sets. 



is illustrated in Table |U We compare the contribution of the full 2n continuum and various 
p-pole parametrizations to the Dirac and Pauli isovector radii of the nucleon [Io| : 

6 f 50M 'lmFy{t) 



(Or=-/ —k^d^ * = 1.2. (8) 

n J AMI 

The first line shows the contribution of the full 2tt continuum from this work to the nucleon 
isovector radius. The given range indicates the error introduced by the different pion form 
factor data sets. The second line shows the p-pole parametrization used in Ref. |2(j. The 
third and fourth lines show the p-pole contribution from Ref. [38| excluding (i) and includ- 
ing (ii) an approximate 2tc continuum, respectively. It is obvious that the simple p-pole 
parametrizations from lines two and three underestimate the full contribution by about 30- 
40% depending on the form factor. The approximate 2n continuum from Ref. 38[ in the 



fourth line comes fairly close for the Dirac form factor but still underestimates the Pauli 
radius by about 20%. 

Inserting the spectral functions into the dispersion relation Eq. (@J), we obtain the 2n 
contribution to the nucleon isovector form factors. The contribution of the spectral function 
in the region t > 50M^ is very small and set to zero in this evaluation. The results for the 
form factors can be fitted by an expression of the form 11, 3 



T// . ai + bj(l — t/ci) 2 I % , s 

fi (t)= qi-t/V) ' I = 1 ' 2 ' (9) 

where a^, 6j, q, and dj are constants. Averaging the results for the three different pion form 
factor data sets [i^, Hl[, the values of the constants are a\ = 1.10788, b\ = 0.109364, 
ci = 0.36963, di = 0.553034, a 2 = 5.724253, b 2 = 1.111128, c 2 = 0.27175, and d 2 = 0.611258. 
The errors in these constants are of the order 4% or less. Using the parametrization from 
Eq. Q, the 2n contribution to the isovector form factors can easily be included in any 
form factor analysis. It is fixed by the pion form factor and ttN data and contains no free 
parameters. 



IV. SUMMARY & CONCLUSION 

In this letter, we have presented a novel analysis of the 2ir contribution to the nucleon 
isovector spectral function using the new high statistics data of the pion form factor by the 
CMD-2 0, KLOE 0, and SND 31] collaborations. The difference in the spectral function 
between the three data sets for the pion form factor is very small. The spectral function 
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displays the contribution of the p peak around t ~ 30M^ and the pronounced enhancement 
close to the two-pion threshold at t — 4M% from the logarithmic singularity on the second 
Riemann sheet. The magnitude of the spectral function increases by about 10% compared 
to the previous analyses. 

The conclusions of Ref. j2?| about the possibility of a long-range pion cloud remain 
unaffected by this change in magnitude: a long-range pion cloud extending as far as 2 fm 
is not compatible with what is known about the 2n contribution to the nucleon isovector 
spectral function. We note that this conclusion might have to be modified if the higher- 
mass pion-continua (3tt, 4tc,...) show a significant threshold enhancement similar to the 
2tt continuum. Given the current state of knowledge, however, this appears unlikely. In 
Ref. |l7J, the threshold behavior of the 3tt continuum was explicitly calculated in heavy 
baryon ChPT and no enhancement was found. Moreover, the inelasticity from four pions in 
7T7T scattering and four-pion production in e + e~ annihilation at low momentum transfer are 
known to be small 34], 3^, H(| . 



Finally, we have calculated the resulting contribution to the nucleon isovector form factors 
and given a convenient parametrization of the result. This contribution is fixed from the 
pion form factor and 7riV-scattering data and contains no free parameters. It can easily be 
used as an independent input in future form factor analyses. This will reduce the number 
of free parameters and ensure the correct spectral function on the nearest part of the cut in 
the time-like region. A new dispersion-theoretical analysis of the nucleon form factors using 



the updated 2ti continuum is in preparation [41 
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